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The charging of a quantum box, coupled to a lead by tunneling through a single resonant level, 
is studied near the degeneracy points of the Coulomb blockade. Combining Wilson's numerical 
renormalization-group method with perturbative scaling approaches, the corresponding low-energy 
Hamiltonian is solved for arbitrary temperatures, gate voltages, tunneling rates, and energies of the 
impurity level. Similar to the case of a weak tunnel barrier, the shape of the charge step is governed 
at low temperatures by the non-Fermi- liquid fixed point of the two-channel Kondo effect. However, 
the associated Kondo temperature Tk is strongly modified. Most notably, Tk is proportional to the 
width of the level if the transmission through the impurity is close to unity at the Fermi energy, and 
is no longer exponentially small in one over the tunneling matrix element. Focusing on a particle- 
hole symmetric level, the two-channel Kondo effect is found to be robust against the inclusion of an 
on-site repulsion on the level. For a large on-site repulsion and a large asymmetry in the tunneling 
rates to box and to the lead, there is a sequence of Kondo effects: first the local magnetic moment 
that forms on the level undergoes single-channel screening, followed by two-channel overscreening 
of the charge fluctuations inside the box. 

PACS numbers: PACS numbers: 73.23.Hk, 72.15.Qm, 73.40.Gk 
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I. INTRODUCTION 

The two-channel Kondo effect 1 ^ is a prototype for non- 
Fermi-liquid behavior in correlated electron systems. It 
occurs when a spin-i local moment is coupled antiferro- 
magnetically to two identical, independent conduction- 
electron channels. Below a characteristic energy scale, 
ksTfc, the system is governed by an intermediate- 
coupling non- Fermi- liquid fixed point, representing the 
fact that neither conduction-electron channel can exactly 
screen the impurity moment. The resulting low-energy 
physics is characterized by anomalous thermodynamic 
and dynamic properties^ Hampering the quest for an 
experimental realization of the two-channel Kondo ef- 
fect is the extreme instability of the non-Fermi-liquid 
fixed point against various perturbations. Any channel 
asymmetry, however small, drives the system to a Fermi- 
liquid fixed point, as does the application of a magnetic 
fields Hence the observation of a fully developed two- 
channel Kondo effect appears hopeless, unless one is able 
to identify a system where the equivalence of the two 
conduction-electron channels is guaranteed by symme- 
try, and all relevant perturbations, such as an applied 
magnetic field, can be tuned to zero. 

One of the leading scenarios for the realization of the 
two-channel Kondo effect is that of a quantum box, either 
a small metallic grain or a large semiconducting quantum 
dot, weakly connected to a lead by a single- mode point 
contact. Near the degeneracy points of the Coulomb- 
blockade staircase, one can map the charge fluctuations 
in the quantum box onto a planner two-channel Kondo 
Hamiltonian, 4 with the two available charge configura- 
tions in the box playing the role of the impurity spin, and 
the physical spin of the conduction electrons acting as a 



passive channel index. The energy difference between the 
two charge configurations corresponds in this mapping to 
an effective magnetic field, which can be tuned to zero 
by varying the gate voltage. Indeed, some signatures of 
the two-channel Kondo effect were recently observed for 
such a setting in semiconductor quantum dots. 5 

However, as recently emphasized by Zarand et al.p 
measurement of the low-temperature, non-Fermi-liquid 
regime of the two-channel Kondo effect sets opposite con- 
straints on the size of the quantum box. On the one 
hand, the charging energy must be sufficiently large in 
order for a measurable Kondo temperature to emerge, 
limiting the box from being too large. On the other 
hand, the mean level spacing in the box must be suf- 
ficiently small compared to fcsTjc, as not to cut off the 
approach to the non- Fermi-liquid fixed point. Hence the 
box cannot be too small. As argued by Zarand et al., 6 
these conflicting limitations cannot be simultaneously re- 
alized in present-day semiconducting devices. The alter- 
nate possibility of using metallic grains is faced with a 
different difficulty of fabricating stable atomic-size con- 
tacts, which are required for obtaining a measurable Tr- A 
Hence the prospects for obtaining a fully developed two- 
channel Kondo effect within Matveev's original picture 
remain unclear. 

In this paper we show that the two-channel Kondo tem- 
perature Tk, and thus the chances for observing the two- 
channel Kondo effect, can be greatly enhanced if tunnel- 
ing between the lead and the box takes place via a single 
resonant level. The study of such resonant tunneling was 
initiated by Gramespacher and Matveev^ who showed 
that one can have a nearly perfect Coulomb staircase, 
even if the transmission coefficient through the impurity 
is one at the Fermi energy. This differs markedly from the 
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case of an energy-independent transmission coefficient, 
where the Coulomb staircase is washed out for perfect 
transmission, & Here we resolve the shape of the Coulomb 
step separating two neighboring charge plateaus, for the 
case of tunneling through a resonant level. 

Using combined analytical and numerical techniques 
we find that the shape of the step is governed at low 
temperatures by the non-Fermi-liquid fixed point of the 
two-channel Kondo effect, similar to the case of a weak 
tunnel barrier^ However, the associated Kondo temper- 
ature is strongly modified. Most notably, Tk is no longer 
exponentially small in one over the tunneling matrix ele- 
ment if the transmission through the impurity is close to 
unity at the Fermi energy, but rather is proportional to 
the width of the level. In general, Tk strongly depends 
on the ratio of the tunneling rates to the box and to the 
lead, which illustrates the inequivalent roles of the two 
rates. If the level is at resonance with the Fermi energy, 
this ratio defines the crossover from weak to strong cou- 
pling. The dependences of Tk on the tunneling rates and 
on the energy of the level are analyzed in detail, as are 
the position and shape of the charge step. 

A potential concern with the above setting has to do 
with the effect of an on-site Coulomb repulsion on the 
impurity level, as it couples the two spin channels. Mod- 
eling the interacting level by a symmetric Anderson im- 
purity, we show that the two-channel Kondo effect is ro- 
bust against the inclusion of a Coulomb repulsion on the 
impurity level, and that Tk is enhanced by a moderately 
large repulsion in the mixed-valent regime. For a large 
on-site repulsion, a local magnetic moment is formed on 
the level. In this regime, Tk decays exponentially with 
on one over the tunneling rates. 

The remainder of the paper is organized as follows: 
Section [H] introduces the physical setting under consid- 
eration. The relation with the two-channel Kondo Hamil- 
tonian is clarified in sec. UTTl for the case of a noninter- 
acting level. An analytic treatment of the weak-coupling 
regime for a noninteracting level is presented in sec. IIVI 
both for a level at resonance and a level off resonance 
with the Fermi energy. This is followed in sec. Ivl bv a 
comprehensive analysis of all coupling regimes using the 
numerical renormalization-group method. The effect of 
an on-site repulsion on the level is studied in turn in 
sec. IVII followed by a discussion and a summary of our 
results in sec. IVIII 
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FIG. 1: Schematic sketch of the physical setting. A quantum 
box and a metallic lead are each coupled by tunneling to an 
intermediate impurity level. The charge inside the quantum 
box is tuned by varying the gate voltage V g , while the energy 
of the level is controlled by varying the gate voltage V%. 



the single-particle dispersion e^B, and by the charging 
energy Ec — e 2 /2C. Here C is the capacitance of the 
box. The mean level spacing inside the box is assumed 
to be considerably smaller than all other energy scales 
in the problem, such that a continuum-limit description 
can be used. The charge inside the box is controlled by 
varying the gate voltage V g , which determines the elec- 
trostatic potential inside the box. We parameterize the 
latter by the dimensionless number N = C g V g /e, where 
C g is the capacitance of the gate, and — e is the electron 
charge. 

Modeling the lead by a noninteracting band with dis- 
persion ekL , the Hamiltonian of the system is given by 

n = J2 e k a cl aa c kaa + e d J2 d Ua (1) 

ol—L.B k,cr <J 
k,tra 

where c\ La {c\ B(J ) creates an electron with spin projec- 
tion a in the lead (box); tkL (ifcs) are the matrix elements 
for tunneling between the impurity and the lead (box); 
and 



k.a 



[ c kBcr 



CkBo 



(2) 



II. COULOMB BLOCKADE WITH RESONANT 
TUNNELING 

The physical setting under consideration is shown 
schematically in Fig. ^ It consists of a metallic lead 
and a quantum box, each coupled by tunneling to an im- 
purity placed in between the lead and box. The impurity 
is assumed to have just a single energy level td in the 
relevant energy range, described by the two creation op- 
erators dX and d\ . The quantum box is characterized by 



measures the number of excess electrons in the box. Here 
and throughout the paper we set the chemical potential 
as our reference energy, namely, all single-particle ener- 
gies are measured relative to the chemical potential. For 
an interacting level, Eq. (JIJ is supplemented by the on- 
site repulsion term 

Hu = Un^fiii, (3) 
where fid a — d\d a are the number operators on the level. 



3 



The Hamiltonian of Eq. features two basic energy 
scales, 



r L = nJ2t 2 kL S(e kL ), 

k 



(4) 
(5) 



corresponding to half the tunneling rates from the impu- 
rity to the lead and to the box, respectively. As shown by 
Gramespacher and Matveevji for Tl^b *C Ec there is 
a nearly perfect Coulomb staircase, even if the transmis- 
sion coefficient through the impurity is one at the Fermi 
energy. The shape of the sharp steps near half-integer 
values of N was left unresolved in Ref. 0, which is the 
objective of the present paper. To this end, we focus 
hereafter on N = n + | + SN, where n is an integer and 
\8N\ <C 1. This range in N corresponds to the step sep- 
arating the two charge plateaus with n and n + 1 excess 
electrons in the box. 

For temperatures well below the charging energy, 
ksT <C Eq, only the n and n+1 charge configurations are 
thermally accessible in the box. Hence one can remove all 
excited charge configurations by projecting the Hamilto- 
nian of Eq. 0} onto the n and n + 1 subspaces. Following 
Matveev^ a spin-i isospin operator S is introduced to la- 
bel the two available charge configurations: S z — 1/2 for 
the n + 1 subspace, and S z — —1/2 for the n subspace. 
The raising and lowering operators, = S x ± iS y , de- 
scribe then transitions between the n and n+1 subspaces, 
corresponding to the addition or removal of a box elec- 
tron. Hence the Hamiltonian of Eq. © is converted to 



where 



U L = E e kLcl La c kLa + e d 

k . a 



J2 d ° d ° 
E*^{ c La+ h - c -} 



(0) 



describes the coupled lead and impurity, 



E 



£kBc\ Bc ,c k B<j ~ eV B S z 



(7) 



with eVe = 2Eq5N describes the isolated box, and 

H tun = E tkB { C IbJ° S+ + H - c -} ( 8 ) 

describes tunneling between the impurity and the box. 
For an interacting level, the above Hamiltonian is sup- 
plemented by the on-site interaction term of Eq. J2J • 

Equations © (JSJ are a straightforward generalization 
of Matveev's original mapping for a weak tunnel barrier: 4 
As in the case of a weak tunnel barrier, the average excess 
charge in the box takes the form 



(9) 



while the capacitance of the junction, 

C(V B ,T) = -d(Q)/dV B , 



(10) 



is proportional to the isospin susceptibility. The connec- 
tion between the Hamiltonian of Eqs. ©~© and the 
two-channel Kondo model is less transparent than for 
a weak tunnel barrier, since the tunneling Hamiltonian 
Htun involves the localized d\ degrees of freedom. As 
we show in the following section, one can still relate the 
Hamiltonian of Eqs. to the two-channel Kondo 

model, by first diagonalizing the quadratic Hamiltonian 
term Ti^ . This gives rise to a new variant of the planner 
two-channel Kondo Hamiltonian, in which the spin-up 
and the spin-down conduction electrons have two distinct 
bandwidths. 



III. RELATION TO THE TWO-CHANNEL 
KONDO HAMILTONIAN 



Much of the underlying physics of the Hamiltonian of 
Eqs. iJfjJl-lJEJ can be understood by converting to a single- 
particle basis that diagonalizes the quadratic Hamilto- 
nian term Hl. The objective of this section is to con- 
struct such a basis. In doing so we assume that the 
level ed lies well within the band (i.e., no bound state 
is formed) , and neglect for simplicity all fc-dependence of 
the tunneling matrix elements t k L and t k s- The latter 
are taken for convenience to be real and positive. 

A convenient representation of the eigen modes of Hl 
involves the c£-electron Green function 



G{z) 



4-2 \ ^ 



1 



LA ^Z~€ kL 



along with the associated phases 

c/) k = arg{G(e fc L - irf)} . 



(11) 



(12) 



Here n is a positive infinitesimal. Introducing the prop- 
erly normalized fermion operators 



^kLa 



"La + tL\G(e kL + irj)\ 

1 



(13) 



4 + *l E 



■f £kL - tk'L + IV 



L fe'Lo 



the Hamiltonian term Hl acquires the diagonal form 

nL = J2 € ^L^kL„ (w) 



fe.fj 



while the d£ operators are expanded as 

4 = ^10(^+^)1^. 



(15) 
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Further converting to the constant-energy-shell operators 

J - 



,t - 

L eBa — 



—==^2S(e-e kL )^l 
J_ ^g(e-e fcfl )4 



(16) 
(17) 



[here Pi(e) and ps(e) are the underlying lead and box 
density of states, respectively] , the full Hamiltonian reads 

H = X! X! / ea \aa a taade - eV B S z (18) 
a=L,B a ■> 

+ lY,J de J de'J(e,e') {at Ba a e , La S+ + K.c.} . 

Here, eVe is equal to 2EcSN; the energy-dependent cou- 
pling J(e, e') is given by 



J(e, e') = 2t L t B y/pB(e) PL (e')\G(e' + irj)\; 



(19) 



and the single-particle operators a ea & obey canonical an- 
ticommutation relations: 



,t 



, j = S(e - e')5 aa '5 a 



(20) 



Note that we have omitted in Eq. I|18|) all those 
conduction-electron channels in both TLl and TLb that 
decouple from the tunneling term TLtun- 

Equation (|18|) should be compared with the cor- 
responding constant-energy-shell representation of the 
planner two-channel Kondo model with a local magnetic 
field. In the latter case, the indices L and B are replaced 
with spin-up and spin-down labels, eVs corresponds to 
the local magnetic field, and a acts as the channel in- 
dex. More significantly, J(e, e') of Eq. i(T§|) is replaced 

with J2Cif(e,e') = y p(e)p(e')J±i where J± is the trans- 
verse Kondo coupling, and p(e) is the joint density of 
states (DOS) of the spin-up and spin-down conduction 
electrons. Thus, identifying the lead and box indices L 
and B with isospin-up and isospin- down labels^ Eq. l(TS|) 
coincides with the planner two-channel Kondo Hamilto- 
nian modulo one crucial difference: the effective DOS for 
the isospin-up and isospin-down conduction electrons are 
markedly different in Eq. I|18|) . While the isospin-down 
DOS is equal to /9_s(e), the effective isospin-up DOS is 
given by the spectral part of G(z)& 

pf{e) = --lm{G(e + lV )} = t 2 L p L (e)\G(e + l r 1 )\ 2 . (21) 

In the wide-band limit, and for e well within the band, 
Eq. 1)2 1|) reduces to the Lorentzian form 



Pf(e) 



(22) 



Hence the effective bandwidth for the isospin-up elec- 
trons is equal to T^, i.e., notably smaller than the isospin- 
down bandwidth Such a large separation of band- 
widths for the isospin-up and isospin-down conduction 



electrons has no analog in the conventional two-channel 
Kondo Hamiltonian, where the two spin orientations are 
identical for a zero magnetic field. Moreover, p^ (e) is 
centered about e d , which corresponds for | | ^ ei- 
ther to a nearly filled band (e d < 0) or to a nearly empty 
band (e d > 0). Below we explore in detail the conse- 
quences of these deviations from the conventional two- 
channel Kondo Hamiltonian, but first let us give some 
heuristic arguments for the expected low-energy physics 
in the case where e d = 0. 

As is well known, the low-energy physics of the planner 
two-channel Kondo model is governed by the dimension- 
less coupling poJ±, where po = p(0) is the conduction- 
electron DOS at the Fermi energy. By analogy with the 
two-channel Kondo Hamiltonian, 



J(0,0) 




(23) 



plays the role of poJ± in the Hamiltonian of Eq. i|18fl . 
Specifically, for a level at resonance with the Fermi en- 
ergy, i.e., e d = 0, Eq. (|53J) is equal to (2/tt) y/T B /T L . 

To the extent that the Hamiltonian of Eq. H18|l still 
flows for e d = and V B = to the intermediate- 
coupling fixed point of the two-channel Kondo effect, dif- 
ferent qualitative behaviors are expected of the Kondo 
temperature Tk in each of the limits Tb <C Tl and 
Tl <C r^. For Tb <C ^l, the bare spin-exchange in- 
teraction i s weak. Hence Tk should be exponentially 
small in \JTlJTb ~ 1/J(0,0). In the opposite limit, 
Tl <C Tb, the dimensionless coupling J(e, e') crosses over 
from weak to strong coupling as |e'| is reduced below 
^bTx ~ \/TjX b- Anticipating a relation between Tk 
and T x , one expects then a Kondo temperature that is 
neither exponentially small in 1 / ^/Tl, nor in 1/ yT b- 

Obviously, the above picture relies heavily on the con- 
jecture that the Hamiltonian of Eq. I|18|) flows for e d = 
and Vb = to the intermediate-coupling fixed point of 
the two-channel Kondo effect, and on intuition borrowed 
from the conventional two-channel Kondo Hamiltonian. 
Although neither assumption is justified a priori, this ten- 
tative picture is shown below to be surprisingly accurate. 



IV. WEAK COUPLING 

We begin our discussion with the limit of weak cou- 
pling, J(0, 0) <C 1, for which an analytical treatment is 
possible. Specifically, we employ a perturbative scaling 
approach based on Anderson's poor-man's scaling^i to 
study two generic cases: (i) e d — and Tb <C Tl, cor- 
responding to an impurity level at resonance with the 
Fermi energy; and (ii) |eci| ^> Tb,Tl, corresponding to 
an impurity level off resonance with the Fermi energy. 
Throughout this paper we assume a symmetric rectan- 
gular form for the underlying lead and box density of 
states, with a single joint bandwidth The latter is 
taken to be much larger than |ed|,r.B, and T^, such that 
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FIG. 2: The basic iterative step in the perturbative scaling 
procedure, for D > Dl- Suppose that the isospin-down (B) 
bandwidth has already been lowered from its initial value D 
to some smaller value D > D' > Dl- The bandwidth D 1 is 
then further reduced to D' — SD' using poor-man's scaling, 
while maintaining the isospin-up (L) bandwidth at Dl- This 
procedure is repeated until D' coincides with Dl , leaving just 
a single common bandwidth for the two isospin orientations. 



p c L s (e) = t 2 L pL(e)\G(e + irj)\ 2 has the Lorentzian form of 



Eq. 123). 



A. Level at resonance with the Fermi energy 

Consider first the case where €d — and Tb <C T l . 
For td = 0, the Lorentzian DOS p c L s (e) is centered about 
the Fermi energy. To proceed with our analytical treat- 
ment, it is convenient to replace pff (e) with a symmetric 
rectangular DOS that preserves both the height of p L S { e ) 
at the Fermi energy, and its total integrated weight: 



(24) 



With this modification, and using the relation J(e, e') 



2ts \I Pb{^)P l S ( £ ')j the Hamiltonian of Eq. I|18|l becomes 



eal La a eL(7 de + ^2 / ^\ B<J a tB ade 

Dl a J-D 

j pD pD L 

+ -yE/_/ e J 1 de'[al Ba a e , L<7 S + +K.c.\ 



eVsSz 



(25) 



where J± = (2 / tt) ^JTb/Tl <C 1 is the effective isospin 
exchange interaction, and Dl is equal to t:Tl/2. Here- 
after we assume that e|Vg| <C D L - 

To cope with the different bandwidths in Eq. 1)25(1 . 
we proceed with perturbative scaling. Using poor-man's 
scaling^i we successively reduce the larger bandwidth 
from D down to Dl, mapping thereby the Hamiltonian 
of Eq. I|25J) onto an effective low-energy Hamiltonian with 
a single joint bandwidth Dl- The basic iterative step 
in this procedure is illustrated in Fig. [3 Suppose that 
the isospin-down (B) bandwidth has already been low- 
ered from its initial value D to some value D' = De~ l , 
< I < \u(D/Dl)- Further reducing the bandwidth to 
D'(l — SI) produces a renormalization to a new interac- 



Dr 



tion term not present in the original Hamiltonian: 

\J2 I de/ de' :al La a e ,La: S z . (26) 

Here : a\ L(J a tl L a ■= a\ L<J a e , La - 6(-e)5(e - e') stands for 
normal ordering with respect to the filled isospin-down 
(L) Fermi sea. Explicitly, A renormalizes according to 
the scaling equation 



dX _ 1 ~ 2 

m - 2 J± ' 



(27) 



where J± = (2/-k)^JT b /T l is the bare isospin-exchange 
coupling in Eq. 125|) . Indeed, the J± interaction term re- 
mains unchanged in Eq. I|25() throughout this procedure, 
apart from the reduced integration range over e. For a 
nonzero e^, there is an additional renormalization of the 
"magnetic" field eVs, discussed below. 

Upon reducing D down to Dl, the new coupling A 
grows from zero to \ J\ kx(D/D L ). Thus, for D' = D L 
one arrives at the effective Hamiltonian 

pD L 

H = / eat aa a eaa de-eV B S z (28) 

a a—L : B Dl 

+ [ D " de [ DL de'{al Ba a t , L aS + + U.c.} 

a J —Dl J — D l 
pD L pD L 

2 a! I de de' \ al LrT a t . L a- a\ B<J a t , B a \ S z 

„ J — Dl J — Dl 



J. 



D L pD, 

de 

-D L J-D 



de I de :al aa .a e > aa : S z , 



where J z = 2V = h J\ \sx{DJDl)- Here we have sepa- 
rated the interaction term of Eq. 1(27(1 into a longitudinal 
isospin-exchange interaction J 2 , and a V term. 

Apart from the extra V term, Eq. 1(28(1 has the form of 
a conventional two-channel Kondo Hamiltonian with an 
anisotropic spin-exchange interaction. As in the case of 
a weak tunnel barrier^ the indices L and B are identi- 
fied in such a mapping with isospin-up and isospin-down 
labels, while the physical spin a serves as a conserved 
channel index. Contrary to the case of a weak tunnel 
barrier, though, the Hamiltonian of Eq. 1281) contains 
a longitudinal Kondo coupling J Zl which can be either 
smaller or larger than J±. For ^J±\n(D/ Dl) <C 1 [i.e., 
(I/^^Tb/Tl ln(D/T L ) < 1] one has J z < J±, while 
for i Jj_ Iii(D/Dl) S> 1 the order is reversed. It should be 
emphasized, however, that Tl (and thus also Tb <C Tl) 
must be exponentially small in this range in order for J z 
to become comparable to J±. 

It is straightforward to verify using perturbative 
renormalization-group (RG) methods that the V term 
acts much in the same way as ordinary potential scatter- 
ing: It is a marginal operator, that does not affect (at 
least not to second order) the Kondo couplings' flow to- 
ward strong coupling. For Vb = 0, the Hamiltonian of 
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Eq. H28|) thus flows to the non-Fermi-liquid fixed point 
of the two-channel Kondo effect, as in the case of a 
weak tunnel barrier^ The corresponding Kondo temper- 
ature can be extracted in turn from known results for 
the anisotropic two-channel Kondo model. In particular, 
for J z <C J± this can be done quite elegantly by iterat- 
ing the standard RG equations backwards (i.e., increas- 
ing the bandwidth D L ), to obtain a planner two-channel 
Kondo Hamiltonian with a bandwidth D* > D L and a 
transverse Kondo coupling J£, that shares the same Tk- 
To leading order in J± and J z / J± one obtains 



D* = \J~Dij5 , Jl = J± 




(29) 



Substituting the above parameters into the expression 
for the Kondo temperature of the planner two-channel 
Kondo Hamiltoniani2i yields 



^bTk = V^bD exp 



7T 

"T 



(30) 



Here, as usual, Tk is given up to a factor of order 
unity, which depends both on the precise definition of 
the Kondo temperature, and on the actual Lorentzian 
form of p c L s (e). 

Two comments should be made about Eqs. I|29 |I -I|H() |I . 
First, as seen in Eq. I|29l) . the effective bandwidth for 
J z <C Jl is neither D nor D L , but rather their geometric 
average, D*. Therefore, the effect of the narrow reso- 
nance that forms on the level is to reduce the effective 
bandwidth in the problem. Second, similar to the case 
of a weak tunnel barrier^ the exponential dependence in 
Eq. (|30|l can be recast in the form exp —ir 2 /2VT , where 

T = 4Tb /Tl is the transmission coefficient through the 
impurity at the Fermi energy for Tb <C Tl. Hence a 
resonant level with = and < Ti acts similar to 
a poorly conducting tunnel barrier with a transmission 
coefficient equal to T — 4Tb/Tl. 



B. Level off resonance with the Fermi energy 

Next we consider a level off resonance with the Fermi 
energy, namely, \e d \ S> Tl,Fb- As noted by Grames- 
pacher and Matveev^ tunneling into and out of the quan- 
tum box are no longer equivalent for ed ^ 0. Depending 
on the sign of ed, this has the effect of either pushing 
down or pulling up the charge plateaus. Using perturba- 
tive scaling we show below that a nonzero ed also shifts 
the position of the degeneracy point, maintaining the 
two-channel Kondo effect at the shifted position of de- 
generacy point. 

To devise a perturbative scaling treatment of the case 
|e<i| 3> r^jTs, one can use either the original Hamil- 
tonian of Eqs. (JBJ— (JSJl, or the equivalent representation 
of Eq. (|18|l . In the former representation, one first re- 
duces the bandwidth from its bare value D to an effec- 
tive bandwidth of the order of le^l, and then performs a 












FIG. 3: Illustration of the density-of-states substitution 
Pl ff ( e ) ~~ * Peff(e), for td < 0. For |e<j| 2> Tl, the Lorentzian 
pff(e) has two prominent features: a narrow peak of width 
~ Tl and weight w ~ 1 centered about td, and a shallow tail 
that crosses the Fermi energy. Each of these features is con- 
veniently mimicked within p e ff(e) by a rectangular structure, 
one narrow and sharp centered about ed, and the other broad 
and shallow centered about the Fermi energy. 



Schrieffer- Wolff transformation 13 to eliminate the charge 
fluctuations on the level. These two steps enter the repre- 
sentation of Eq. i|18fl hi a unified fashion through the en- 
ergy dependence of the coupling J(e, e'), which is sharply 
peaked as a function of e' at ed- 

We have carried out both the perturbative scaling 
approach based on the Hamiltonian representation of 
Eqs. 10-©, and the scheme based on the Hamiltonian 
representation of Eq. I|18fl . Both procedures give the 
same results to leading order in r j, / 1 | and r# / 1 £d | (note 
that the Schrieffer- Wolff transformation is designed to 
capture only the leading order in these parameters) . For 
the sake of consistency with the analysis of the previous 
subsection, we present below the approach based on the 
Hamiltonian representation of Eq. (|18fl . 

Similar to the case of a level at resonance with the 
Fermi energy, it is convenient to replace p c L s (e) with a 
simplified density of states that captures the essential 
features of p L s (e), and allows for an analytical treat- 
ment of the problem. Specifically, p L S (e) has two promi- 
nent features: a narrow resonance of width ~ T L and 
weight 1 — O (r^/ledl) centered about e = ed, and a shal- 
low tail that crosses the Fermi energy and provides the 
low-energy excitations for the development of the Kondo 
effect. To mimic these two features, we replace p c L s (e) 
with the double-rectangular density of states illustrated 
in Fig. 



Peff(e) 



w 
2D~l 



6(D L -\e-e d \) + ^6{D ri 
ne^ 



|e|). (31) 



Here D m = \ed\ — D L is a crude measure of the extent 
of the tail that crosses the Fermi energy, D L — ttTl/2 
corresponds to half the width of the resonance at e = e.d, 
and w = 1 — 2Tj J D m /(Tre d r ) « 1 is the effective weight of 
the resonance. 

Obviously, there is some arbitrariness in our choice of 
Peff(e), which differs in details from p c L s (e). In particu- 
lar, the extent of the tail that crosses the Fermi energy 
is somewhat exaggerated, while the opposite tail (the 
one extending away from the Fermi energy) is absent. 
Nevertheless, this choice of p c s (e) is compatible with the 
approximations made in the Schrieffer- Wolff transforma- 
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tion, and gives the correct results to leading order in 
T L /\e d \ and T B /\e d \. Substituting p eff (e) in for p e L s (e), 
the Hamiltonian of Eq. Ijl8[) becomes 



U = I £al La a eLa de- eV B S z 
+ / ea lB* a tBade 

„ J-D 

f-D pD p 

V / de de'V^ 7 ) 
; J-d J-d p 

x { a I B(T a €'L<TS ,+ + H.c. j 



(32) 



where D p = \e d \ + DL- Hereafter we assume that e|Ve| <C 

To treat the Hamiltonian of Eq. (|32[) . we proceed with 
perturbative scaling. This is done in two stages. First 
the larger bandwidth is successively reduced from its bare 
value D down to D p , leaving just a single common band- 
width for the two conduction seas. This common band- 
width is subsequently reduced from D p down to D m . 

The first step in the above scheme is similar to the one 
implemented in the previous subsection, when treating 
a level at resonance with the Fermi energy. There are, 
however, three important modifications. First, the new 
interaction term generated upon scaling has the form 



D„ pD 



A } ' I de de'^p cS (e)p cS (e') : a\ Liy a c > Lr7 : S z , 

a J - D P J - D P 

(33) 

which differs from Eq. I|2t)l) in the extra square roots of 
p e f{ that enter the integrand. Second, the coupling A 
renormalizes according to the scaling equation 



dX 2T r 



dl 7T ' 



(34) 



which likewise differs from Eq. (|27|) . Lastly, the voltage 
Vb is renormalizcd at T = according to 



dV B . 2r, 

— = sign(e d ) 

dl 7T 



l + e l \e d \/D 



(35) 



D 



Here we have distinguished the running parameter Vb 
from its bare value Vb, and omitted higher order correc- 
tions in Fi/|ed|. Upon reducing the larger bandwidth 
from D down to D p , the coupling A thus grows from 
zero to (2T b / tt) ln(D / Dp) , while eV B evolves from eVs 
to eV B + sign(e d )(2r B /tt) In 

Since we are interested in |e<j| ^> Tl, one can pro- 
ceed to eliminate all excitations in the energy range 
D p > e > D m in one step, by working with a finite 
51 = 2Dl/D p <C f. Within the Hamiltonian represen- 
tation of Eqs. this step is equivalent to carrying 
out the Schrieffer- Wolff transformation. At the end of 
this procedure one arrives at an effective Hamiltonian of 



the form of Eq. J33J), where Dl is replaced with D m , plus 
two additional potential-scattering terms: 

pD m *D m 

V ±J2j de J de ' {■ a lLa a ^'La- ± ■ a\ B<j a t , B a ■] ■ 

(36) 

The effective coupling constants entering the resulting 
Hamiltonian are quite different, however, from those in 
Eq. I|28|l. and are given by 



- 2 vrzr^ 

J± i — i — ) 

j 

• J z — I I j 

2tt |e d | ' 

v+ = -^u slgn{ed) ' 

= 47R slgn(ed) ' 



(37) 
(38) 
(39) 
(40) 
(41) 



Here we have omitted corrections that are smaller by fac- 
tors of r a /|ed| or (T a /n\e d \) \n(D/\e d \) than the leading- 
order terms. In addition, the voltage V B is renormalized 
at T = according to 



2r / d 

eV B -> eV B = eV B + sign(e d ) — -In ( f + - — - 



(42) 



It is straightforward to verify using either poor-man's 
scaling or bosonization (in combination with a canon- 
ical transformation) that the V± terms are marginal, 
and do not affect the zero-temperature fixed point of 
the Hamiltonian other than through an additional shift 
of Vb- Specifically, neglecting the renormalization of 
J z , Eq. (|42|l acquires the additional small correction 

eV B -> eV B + D p 161n(2) (vV+ + J z V-/2\. Hence the 

system continues to undergo the two-channel Kondo ef- 
fect for | erf | 3> Tl,F b , albeit at a shifted position of the 
degeneracy point, approximately given by 



eV 2 cK = 



-sign(e d ) In 



D 

M 



1 + . (43) 



It should be noted, however, that the associated Kondo 
temperature is quite sensitive to the ratio of Tl to Tb, 
which fixes the ratio of J± to J z in Eqs. l(3*7 )l -[|3*% )l . For 
Tb, one recovers the exponential form Tk oc 



exp 



-7T 2 /2\/T 



where T = 4TLT B /\e d f is the trans- 
mission coefficient through the impurity at the Fermi en- 
ergy. By contrast, the Kondo temperature depends in a 
power-law fashion on fx, for Tb S> Tl ■ 



Th V 'ft 



7r|e d |/2r f 



(44) 



Regardless of the ratio Tl/Tb, Tk decays exponentially 
with \e d \, for \e d \ > r L ,r B . 
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V. GENERAL COUPLING 

Based on perturbative scaling, our treatment thus far 
was confined to the weak-coupling regime, J(0, 0) <C f. 
We now turn to a nonperturbative study of all parame- 
ter regimes, ranging from weak to strong coupling. To 
this end we go back to the Hamiltonian of Eqs. 
and employ Wilson's numerical renormalization-group 
(NRG) method. 14 Originally developed for treating the 
single-channel Kondo Hamiltonian^ this nonperturba- 
tive approach was successfully extended to the Ander- 
son impurity model (both the symmetric*^ and asymmet- 
rioi& models), the two-channel Kondo HamiltonianjiLiSi 
different two-impurity clustersji2i22i2i and a host of re- 
lated zero-dimensional problems. Below we adapt this 
approach to the Hamiltonian of Eqs. (JSJ-lJSJ. 



A. The Numerical renormalization group 

At the heart of the NRG approach is a logarithmic 
energy discretization of the conduction band around the 
Fermi energy. The conduction electrons within each en- 
ergy interval [— DA~' l , — DA - ™ -1 ] and [DA - ™ -1 , DA - ™], 
n = 0, 1, 2, • • •, are replaced by a single degree of free- 
dom per spin orientation, such that all energy intervals 
contribute equally to the infra-red divergences which are 
immanent in the problem. Here A > 1 is a discretiza- 
tion parameter, with the full Hamiltonian recovered for 
A — * 1 + . Using an appropriate unitary transformation, 14 
the conduction band is mapped onto a semi-infinite chain, 
with the impurity coupled to the open end, and the hop- 
ping matrix elements decreasing exponentially along the 
chain. For the problem at hand there are two separate 
bands, one for the lead and one for the box. Hence four 
different Wilson shell operators are required at each point 
along the chain: /^ CTT j, where a = L,B labels the band 
(lead or box), a =f, J, is the spin index, and n = 0, 1, 2, • • • 
enumerates the position along the chain. In this manner, 
the full Hamiltonian of Eqs. is recast as a double 

limit of a sequence of dimensionless NRG Hamiltonians: 



n 



lim lim 

A^l+ JV-»a 



{d a a- 



(N-l)/2 n 



A 



with D A equal to D(l + A) -1 /2, and 



n 



A 



JV-l 
A 2 



D A 



u 

Dl 



:n d y. :n dl : 



(45) 



(46) 



E {wL<As+ + hflaod* + h.c} 



Physically, the shell operators /2 ct0 represent the lo- 
calized states in each band, to which the impurity level 
is directly coupled. The subsequent shell operators f^ an 
correspond to wave packets whose spatial extent about 
the level increases approximately as A™/ 2 . All informa- 
tion on the underlying band structure is contained in the 
hopping coefficients $„ a p which are obtained from ap- 
propriate integrals of the density of states^ In this paper 
we use a symmetric rectangular density of states for both 
the lead and the quantum boxjSi for which one has the 
explicit expression 1 ^ 



1 - A-(" +1 ) 



y/(l ~ A-( 2 "+ 1 ))(l - A-( 2 ™+3)) 



(47) 



Although the discretized form of the Hamiltonian is 
exact only in the limit A — > 1 + , in practice one works 
with a fixed value of A > 1. As shown by Wilson, 14 
the error introduced by not implementing the limit A — > 
1 + is perturbative and small. Starting with the local 
Hamiltonian 7io, the sequence of NRG Hamiltonians are 
iteratively diagonalized using the NRG transformation 

'H N+1 =VXH N + J2^No l {fL N+1 f a aN+'R.C.y (48) 



At each iteration, four new shell operators are intro- 
duced, enlarging the Hilbert space by a factor of 2 4 = 16. 
Since it is numerically impossible to keep track of such an 
exponential increase in the number of basis states, only 
the lowest N s eigenstates of 7Ya are retained at each iter- 
ation. These N s states are used in turn to construct the 
eigenstates of Hn+i using Eq. (|48|l . Thus, two approxi- 
mations are involved in the NRG algorithm: discretiza- 
tion and truncation. Each of these approximations can 
be systematically controlled by varying A and N s . 

The eigenstates of Hat so obtained, Hn\1)n = Ef\l) a, 
are expected to faithfully describe the spectrum of the 
full Hamiltonian H on the scale of D N = DaA"^ 1 )/ 2 . 
Hence, they can be used to compute thermodynamic av- 
erages at the temperature A^Ta = Da//3, where (3 is a 
numerical factor of order unity^ Specifically, the ther- 
modynamic average of an observable O at the tempera- 
ture T/v is approximated by 



(6) 



N 



E 



-0K 



J N 



n(1\0\1) n , 



where 



Z 



N 



(49) 



(50) 



A-l 



E E A_t ^ {/Ln+l/«m + H.C.} 



n— aa 



Here the t a 's are related to the hybridiz ation widths of 
Eqs. (01)-© through D A t a = ^2T a Dj-K, while the pref- 
actor A^ -1 " 2 guarantees that the low-lying excitations 
oi Hn are of order one for all N. 



In this way, one can compute thermodynamic averages 
at a decreasing sequence of temperatures. Note that the 
effect of truncation in Eqs. I|49 p — l|5U| ) can be systemat- 
ically reduced by increasing the number of NRG states 
retained at each iteration. 

Apart from the local Hamiltonian 7i , which involves 
the extra d a degrees of freedom, the NRG formulation of 
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eV B /k B T K 



FIG. 4: Smearing of the charge step at different temperatures. 
Here 6Q = (Q) + e(n + ^ ) is the excess charge inside the box, 
measured relative to the mid point between the two charge 
plateaus, fcsTif / D — 0.0014 is the corresponding two-channel 
Kondo temperature, and Tl/D — Vb/D equals 0.0039. The 
number of NRG states retained is equal to 2300, and A = 2.3. 
At temperature T, the charge step is smeared over a range of 
eVs ~ max{fcsT, fesTif}. The slope of the step at Vb — 
continues to steepen with decreasing temperature (see inset), 
diverging logarithmically as T — > 0. 



our problem is equivalent to that of the anisotropic two- 
channel Kondo Hamiltonian>i& It remains so also for an 
interacting level. Similar to the anisotropic two-channel 
Kondo Hamiltonian, the Hamiltonian of Eqs. ©-JSJ pos- 
sesses three underlying symmetries: SU (2) channel sym- 
metry (spin symmetry), conservation of the total elec- 
tronic charge, and conservation of the z component of 
the total isospin operator: 



Si"*" 1 = \ £ (c{ La C kLa - c{ Ba C kBa ) + l - £ 4da+S x . 

(51) 

Each of the NRG Hamiltonians Hat is block-diagonal in 
the conserved quantum numbers, enabling a reduction in 
the size of the matrices to be diagonalized. In our code 
we exploited only the conservation of the total electronic 
charge and the z components of the total isospin and 
physical spin, ignoring the full SU(2) spin symmetry of 
the problem. This necessitated keeping a larger number 
of NRG states, typically around N s = 2300. 

Note that it is straightforward to include a finite on- 
site repulsion U within the NRG, as it only enters the 
local Hamiltonian H.q. Apart form modifying the eigen- 
energies and eigenstates of Ho, a finite Coulomb repul- 
sion U has no effect on the formulation of the NRG. 
Hence treatment of an interacting level is computation- 
ally equivalent to that of a noninteracting level. 



FIG. 5: The capacitance C(V B ,T) versus T, for Vl/D = 
Vb/D = 0.0157 and different values of Vb. Here k B T K /D = 
0.0063, A = 2.3, and the number of NRG states retained 
is equal to 2300. The deviation in voltage from the degen- 
eracy point, eVe = 2EcSN, takes the values eVs /IsbTk = 
0,0.01,0.03,0.1,0.3, and I, according to the labels attached 
above each individual curve. 

B. Level at resonance with the Fermi energy 

1. Two-channel Kondo effect 

We begin with a noninteracting level at resonance with 
the Fermi energy, i.e., e<j = U = 0. Due to particle- hole 
symmetry, the location of the degeneracy point is pinned 
in this case at V B = 0. 

Figure 0] shows the temperature dependence of the 
charge step, for T L /D = Y B /D = 0.0039. Here and 
throughout the paper we parameterize the charge step 
by the excess charge inside the box, measured relative to 
the mid point between the two charge plateaus: 

SQ(V B ,T) = {Q)+e(n+^\. (52) 

For ed — and equal hybridization widths, T L = T B , the 
parameter J(0, 0) in Eq. I|23() equals 2/tt m 0.64, which 
lies beyond the perturbative scaling analysis of sec. IIVI 

As seen in Fig. 0J the charge step at temperature T is 
smeared over a range of eV B ~ max{fcsT, k B Tj(}- Here 
ksTx/D = 0.0014 is a new low-energy scale, the Kondo 
temperature, whose precise definition is given below. In 
accordance with Matveev's scenario for a weak tunnel 
barrier^ the slope at V B = continues to steepen with 
decreasing T, consistent with the development of a two- 
channel Kondo effect. Indeed, the line shapes in Fig.0]are 
quite similar to those obtained for a weak tunnel barrier 
using the noncrossing approximation^ 

The emerges of the two-channel Kondo effect for such 
intermediate coupling is confirmed by the capacitance 
line shapes, C(V Bl T) — —d(Q)/dV B , which are plot- 
ted in Fig. for T L /D — T B /D = 0.0157^i At the de- 
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FIG. 6: Zero-temperature smearing of the charge step. Open 
circles are the calculated NRG points; the full line shows the 
analytic formula of Eq. 1551 with / = 1.85. The latter ex- 
pression for SQ(Vb,0) relies on the logarithmic fit C(0, T) = 
{e 2 /20k B T K )[\n(T K /T) + B], with k B T K /D = 0.0014 and 
B = 2.16. The quality of this fit for C(0, T) is demonstrated 
in the inset. Here open circles are the NRG data points, while 
the full line shows the logarithmic fit. All model and NRG 
parameters are the same as in Fig. 0] 



generacy point, Vb = 0, the capacitance diverges loga- 
rithmically with decreasing temperature, in accordance 
with the characteristic ln(T) divergence of the mag- 
netic susceptibility in the two-channel Kondo effect £ For 
e|Vs| <C UbTk, this logarithmic temperature dependence 
is cut off about an order of magnitude below the isospin 
polarization scale 



T sp = {eV B f/k 2 B T K , 



(53) 



which governs the crossover from non-Fermi-liquid to 
Fermi-liquid behavior (Pauli-like susceptibility). Thus, 
the effect of eVg is identical to that of a magnetic field in 
the two-channel Kondo effect^ Indeed, the capacitance 
line shapes of Fig. [5] are very similar to the susceptibil- 
ity curves obtained by the Bethe ansatz for the isotropic 
two-channel Kondo model^ with eVs playing the role 
of an applied magnetic field. 

The close resemblance with the Bethe ansatz curves 
for the magnetic susceptibility of the two-channel Kondo 
model provides us with a precise procedure for defin- 
ing the two-channel Kondo temperature Tk- Specifically, 
throughout this paper we define the two-channel Kondo 
temperature by the Bethe ansatz expression for the slope 
of the ln(T) diverging term in the zero-field susceptibil- 
ity^ which translates in this case to 



C(0,T) 



20k B T K 



HTk/T). 



(54) 



- (2r B /D)"exp[- Jt (r L /r B )""/4] 

o NRG: N s =1400 
■ NRG: N s =1500 
A NRG: N =2300 
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FIG. 7: The Kondo temperature Tk versus Tl/^b, for 
T L /D = 0.098 and T L > r fl . Here A = 2.3, while the 
number of NRG states retained is varied from N s = 1400 
to N a = 2300. For F_g <C Tl, the Kondo temperature is well 
described by the weak-coupling formula of Eq. I|30[l . with a 
pre-exponential factor equal to y/2 (solid line) . 



are fitting parameters. The quality of our fits is demon- 
strated in the inset to Fig. [U for the same set of model 
parameters as used in Fig. ^ 



2. Shape of the charge step 

As is evident from the curves of Fig. \5\ the satu- 
rated zero-temperature capacitance for < eVs < ksTx 
closely follows the zero-voltage capacitance C(0, T) at 
the crossover temperature T sp . This suggests that one 
can approximate C(Vb, 0) at small voltages by C(0, T = 
fT sp ), where T sp is defined in Eq. lj5"3*j) . and / is a nonuni- 
versal constant of order unity. Using our logarithmic fit 
for C(0,T), one can then integrate C(0,T = fT sp ) with 
respect to Vb , to obtain the following analytic expression 
for the shape of the zero-temperature charge step: 



6Q(V B ,0) 



e 2 V B 
20k B T K 



\k B T K 



(55) 



Thus, to extract Tk we fitted the capacitance C(0, T) to 
the form (e 2 /20k B T K ) \\n(T K /T) + B], where T K and B 



with B' = B + 2-\n{f). 

Figure El shows a comparison of Eq. I|55|l with the zero- 
temperature charge step obtained from the NRG, for the 
same set of model parameters as used in Fig. 0J The 
NRG data points were obtained by going to a sufficiently 
low temperature, such that 6Q has saturated for each 
voltage point displayed at its T — > value. As seen in 
Fig. El Eq. with / = 1.85 well describes the shape of 
the zero-temperature charge step up to e|Vs| ~ ksTx- 
For e\V B \ > ksT^ there are visible deviations from the 
NRG data points, which stem from the breakdown of the 
logarithmic approximation for C(Vb,0). 
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FIG. 8: The Kondo temperature Tk versus \/Fb, for Fl/D = 
0.0039. Here A = 2.3, while the number of NRG states re- 
tained is equal to 2300. The Kondo temperature monoton- 
ically increases with r_g, initially rapidly for Tb < Fl, and 
then more moderately for Fb > Fl- The value Fb = Fl 
is marked by an arrow. Inset: Convergence of Tk with the 
number of NRG states retained N s , for Fb = 4Tl. Below 
N„ — 2000, finite-size effects are relatively large at this point. 
Above N s = 2100, the Kondo temperature is converged to 
within 1%. 



3. Low-energy scale 

So far, we have focused our attention on intermedi- 
ate coupling, Tr, — Tb- Varying Tl/Tb, we confirmed 
that the two-channel Kondo effect persists for all ratios 
of Tl/Tb, ranging from weak (Tb <C Tl) to strong 
(r^ ^> Tl) coupling. The main effect of Tl/Tb is to 
modify the Kondo temperature Tk, as described below. 

Figure[7|shows Tk versus Tl/Tb in the weak-coupling 
regime, Tb < Tl- Here Tl/D is kept fixed at 0.098, while 
is reduced as to increase Tl/Tb- The Kondo tem- 
perature was extracted from the Bethe ansatz expression 
of Eq. (|54fl . For Tb <C Fl, excellent agreement is ob- 
tained^ with the perturbative scaling result of Eq. (f 301) , 
up to a pre-exponential factor equal to \J~2. This con- 
firms the exponential dependence of Tk on ^/Tl/Tb- 
Deviations from Eq. l|30[l are observed as Tb approaches 
Tl, signaling the crossover to intermediate coupling, and 
the breakdown of perturbative scaling. Note that the 
NRG results are practically free of finite-size effects; only 
small variations are seen upon going from N s — 1400 to 
N s = 2300. 

An extension of Fig.Oto intermediate and strong cou- 
pling is presented in Fig. [S] Here a smaller value of 
Tl/D = 0.0039 was chosen, as to enable larger ratios 
of Tb to Tl, while maintaining Tb *C D. Fixing Tl, the 
Kondo temperature monotonically increases as a function 
of r#, revealing two distinct regimes. For weak coupling, 
Tb <C Fl, one recovers the exponential dependence of 




| (2r L /7i.D) 1/2 

FIG. 9: The Kondo temperature Tk versus \ZTl, for Fb/D = 
0.0039. Here A = 2.3, while the number of NRG states re- 
tained is equal to 2300. Contrary to Fig. |H] where Tl is kept 
fixed, the Kondo temperature increases with increasing Fl 
for r_L < Fb, but then decays exponentially with \/Fl for 
Fl > Fb- In between, Tk has a maximum for Fl ~ Fb (the 
point Fb = r_L is marked by an arrow). Inset: A plot of Tk 
versus Fl, for Fl < Fb- 

Eq. (|30(l . For Tb 3> Tl, this exponential form crosses 
over to an approximate linear dependence on y/Ts- In 
the latter regime, Tk is no longer exponentially small in 
one over the tunneling matrix elements, as we expand 
below. 

In Fig. E| the ratio Tl/D = 0.0039 is kept fixed and 
is varied. Figure |5l shows the complementary picture, 
whereby Tb/D = 0.0039 is kept fixed and Tl is varied. 
As is clearly seen from comparison of Figs. 00 and EJ the 
two hybridization widths Tl and Tb have inequivalcnt 
roles in the two-channel Kondo effect that develops. In 
particular, Tk increases monotonically as a function of 
Tb for fixed Tl, but depends nonmonotonically on Tl 
when is kept fixed. In the latter case, Tk initially 
grows with Tl for Tl < Tb, but then decays exponen- 
tially with ^/Tl for 1^ > Tg. In between, Tk has a 
characteristic peak for Tl ~ Tb- 

Experimentally, one can test the predictions of Figs. [HI 
and by separately tuning Tl and Tb using appropriate 
gate voltages. Perhaps the most natural parameteriza- 
tion of the combined lead-level-box junction, though, is 
in terms of the transmission coefficient through the im- 
purity at the Fermi energy: T = 4r L r B /(r L +r B ) 2 . 
Here we have set td = in the expression for T, which 
depends symmetrically on the hybridization widths Tl 
and r^. Since Tl and Tb have inequivalcnt roles in de- 
termining Tk, it is clear that the Kondo temperature 
depends differently on T, for Tl > Tb and Tl < Tb- 

In Fig. we rescaled the Kondo temperatures of 
Figs. |H1 and by the level broadening Tl+Tb, and plot- 
ted the resulting ratio as a function of T. Even after 
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FIG. 10: The Kondo temperatures of Fig. [S] (where Yl is 
kept fixed) and Fig. (where Tb is kept fixed), rescaled by 
the level broadening Tl + Tb, and replotted as a function of 
the transmission through the impurity at the Fermi energy, 
T = 4ri,r s /(r L +r s ) 2 . The left-hand side of the graph cor- 
responds to Tb < Tl, while the right-hand side corresponds 
to T L < T B - As a function of T, the ratio k B T K /(r L + Tb) 
has a peak near perfect transmission. For T « 1, the Kondo 
scale roughly equals one fifth of the level broadening. 



separating the regimes Tl > Tb and < Tb, the 
ratio UbTk/^l + Tb) is not an exclusive function of 
T. Rather, it depends on the way in which T is tuned. 
Nevertheless, the resulting T dependence has several dis- 
tinct characteristics. Near perfect transmission, the ratio 
^bTk/^l + Tb) is greatly enhanced. For Tb < Tl, this 
ratio is peaked at T — 1, while for T L < Tb the peak 
is slightly shifted below T = 1. Most notably, UbTk 
roughly equals §.2(Tl + T b) near perfect transmission, 
which is many-fold larger than the Kondo temperature 
obtained for a tunnel barrier with comparable tunneling 
matrix elements. 

Indeed, while Tk is exponentially small in one over 
the tunneling matrix element in the case of a weak tun- 
nel barrier^ for Tl = Tb it depends approximately lin- 
early on the level broadening. This important point is 
demonstrated in Fig. ^2 where Tk and HbTr /(J'l + ^'b) 
are plotted as a function of Tl — Tb- For over two 
decades in (r^ + Tb), the Kondo scale UbTk roughly 
equals Q.2{Tl + Tb), establishing the departure from the 
familiar exponential form of Tk- This dramatic enhance- 
ment of Tk should have important implications for the 
observation of the two-channel Kondo effect in semicon- 
ducting devices. Tuning rj, w 1^ to be smaller than the 
charging energy yet notably larger than the level spacing 
inside the box, one can possibly realize a Kondo temper- 
ature that significantly exceeds the level spacing, thereby 
opening the door to the observation of the two-channel 
Kondo effect in semiconducting devices. 

It should be noted that finite-size effects are small in 
Fig. El down to (T L +T B )/D ~ 5 x 1CT 4 (compare 
N s = 2300 with N s = 2500), but rapidly increase be- 
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FIG. 11: k B T K versus T L + T B , for T L = T B (i.e., perfect 
transmission through the impurity at the Fermi energy). For 
over two decades in (Tl + Tb)/D, the Kondo scale UbTk 
depends approximately linearly on Tl + Tb, with a propor- 
tionality factor roughly equal to 0.2 (see inset). 



low this value (not shown). Therefore, we are unable to 
reliably determine whether or not the approximate linear 
dependence of Fig. 1111 persists down to smaller values of 
(T l +Tb)/D. 



C. Level off resonance with the Fermi energy 



As discussed in sec. IIVBI the position of the degen- 
eracy point is no longer pinned at Vb = for nonzero 
Cd- Specifically, for €d ^ and ksT < |e<j|, the aver- 
age occupation of the level deviates from one-half per 
spin orientation. Depending on the sign of 6d, the level 
is more available then cither for tunneling into (e^ < 0) 
or out of (e e i > 0) the quantum box, which generates a 
dynamic "magnetic" field acting on the isospin. At the 
degeneracy point, this dynamic field must be balanced 
by a nonzero Vb, causing a shift in the position of the 
degeneracy point. The resulting low-energy physics at 
the shifted degeneracy point is that of the two-channel 
Kondo effect, with a reduced Kondo temperature that 
decays exponentially with \ed\, for Tl,Tb <C |ej|. 

The above picture was established in sec. IIVBI using 
perturbative scaling. We confirmed its validity using the 
NRG. In particular, Fig. 1121 shows the zero-temperature 
charge step [i.e., SQ(Vb,0) versus Vb], for Tl = Tb and 
different values of > 0. The corresponding charge steps 
for ed — * — £d are obtained by simultaneously reflecting 
the curves of Fig. H2I about V B = and SQ = 0M As 
before, the NRG curves were computed by going to a 
sufficiently low temperature, such that SQ has saturated 
at its T — > value for each voltage point displayed. 

With increasing td > 0, the charge step initially shifts 
toward more negative voltages, before moving back in di- 
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FIG. 12: Evolution of the charge step with increasing e^, for 
F L /D = T B /D = 0.01. Here A = 2.3, while the number of 
NRG states retained is equal to 2300. With increasing > 0, 
the charge step initially shifts toward more negative voltages, 
before moving back in direction of Vb = 0. As exceeds Tl, 
there is a significant narrowing of the step separating the two 
adjacent charge plateaus. For e<j — Tl, the position of the 
degeneracy point V2CK is indicated by an arrow. 

rection of Vb = 0. The width of the step separating the 
two adjacent charge plateaus becomes narrower with in- 
creasing ed, in accordance with the notion of a decreasing 
Tk- This narrowing of the step is particularly transpar- 
ent for larger values of €d (exemplified by edjT L = 10 
in Fig. ll2f> . when the underlying Kondo temperature be- 
comes exponentially small in ed/T L . Also apparent is the 
development of a slight asymmetry in the shape of the 
step. As a function of Vb, the approach to the left charge 
plateau with n excess electrons in the box is more rapid 
than the approach to the right charge plateau with n + 1 
excess electrons. 

One can understand this asymmetry in the shape of the 
charge step along the following lines. Consider the limit 
(d > (Tl + Th), such that the level is essentially empty. 
For values of Vb sufficiently removed to the right of the 
step, the box is predominantly in the n + 1 charge con- 
figuration. The weight of the n charge configuration can 
then be estimated by simple perturbation theory with re- 
spect to ts , which yields a contribution of order t B . Here 
we made use of the fact that the empty level is free to 
absorb a box electron. By contrast, for values of Vg suf- 
ficiently removed to the left of the step, the empty level 
has no electron to donate to the box, which blocks any 
direct matrix element between the n and n+1 charge con- 
figurations. Hence the level must first be charged with 
a lead electron before this electron can be passed on to 
box. This reduces the weight of the n+1 charge con- 
figuration by an extra factor of order t\, rendering the 
n charge configuration more stable. Therefore, charge 
fluctuations are more prominent to the right of the step. 

A similar effect, albeit more pronounced due to the 




FIG. 13: The capacitance at the shifted degeneracy point, 
for <L d = Y L = Ts = 0.01£>. Here k B T K /D = 0.0042, 
A = 2.3, and N s — 2300. As indicated by the arrow in 
Fig. 1121 the degeneracy point is slightly shifted away from 
the position of the mid-charge point (cV2Ck/Tl = —1.4 ver- 
sus eV* /Tl = —1.22). For V = V2CK, the capacitance 
diverges logarithmically with T — > 0, in accordance with 
the onset of the two-channel Kondo effect. This is demon- 
strated by the solid line, which shows the logarithmic fit 
C(V 2C k,T) = (e 2 /20k B T K ) [ln(T K /T) + B] with k B T K /D = 
0.0042 and B = 2.48. Inset: Temperature dependence of 
the charge curve near the degeneracy point. The solid, dot- 
ted, dotted-dashed, dashed, and thin-solid lines correspond to 
T/T K = 4xl0" 9 , 0.042, 0.15, 0.34, and 0.78, respectively. The 
position of eVzcK is indicated by an arrow. 



larger values of T used, is seen in the second-order per- 
turbation theory of Gramespacher and Matveev^ who 
found that the charge plateaus are "pushed down" for 
td > 0. Due to the breakdown of perturbation theory, 
these authors were unable to access the vicinity of the 
degeneracy point. Here we can exploit the NRG to quan- 
titatively trace the shift in the position of the charge step. 
Explicitly, we focus below on the mid-charge point V* , 
defined as the voltage for which 5Q(V* , 0) = 0. 

An important aspect of the asymmetric line shape for 
6d 7^ is a separation between the mid-charge point V* 
and the degeneracy point V2CK, defined as the voltage 
at which two-channel Kondo physics emerges. In the 
presence of particle-hole symmetry, V* and V2CK coin- 
cide. This is no longer the case away from particle- 
hole symmetry, as demonstrated in Figs. ^| and ^| for 
6d = T L = Tb = 0.01Z3. Here the degeneracy point is 
shifted further away from Vb = than the mid-charge 
point (eV 2C K/T L = -1-4 versus eV*/T L = -1.22). 
While the mid-charge point is detected from the con- 
dition SQ(V*,0) = 0, the degeneracy point is identified 
as the point where the low-temperature charge curves in- 
tersect (see inset to Fig.^Jl. As seen in Fig. the ca- 
pacitance diverges logarithmically with decreasing tem- 
perature at Vb = V2CK , confirming the onset of the two- 
channel Kondo effect. While clearly not that of a lo- 
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cal Fermi liquid, the finite-size spectrum at Vb = V 2 ck 
deviates from that of the standard two-channel Kondo 
Hamiltonian, which is likely due to the way in which 
Ed 7^ breaks particle-hole symmetry in the Hamiltonian 
of Eqs. Although V* and V 2 qk differ for nonzero 

ed, they do closely follow one another. We therefore con- 
centrate hereafter on the mid-charge point, which is much 
easier to trace numerically. 

Figure El depicts V* as a function of td > 0, for 
T L /D = Tb/D — 0.01. As indicated in Fig. [TJ V* 
has a minimum at an intermediate energy 2.5r^,, 
reaching a minimal value of roughly —2T^. For 3> T^, 
V* gradually increases to zero according to the analytic 
formula of Eq. I|43|l. The latter expression (which was de- 
rived, strictly speaking, for the degeneracy point V 2 cyl) 
is shown for comparison by the solid line in Fig. 1141 For 
£(2 ^> Tz, , there is good agreement between the NRG and 
perturbative scaling. For Ed ^> -D, the position of the 
mid-charge point (as well as that of the degeneracy point) 
once again approaches V* — 0, due to the suppression of 
charge fluctuations on the level. 

Two comments should be made about the predictions 
of Figs. El an d El First, one can experimentally test 
these predictions by tuning the gate voltage Vi in the set- 
ting of Fig.^ which has the effect of varying td- Second, 
these predictions are for zero temperature. In general, 
the position of the mid-charge point is temperature de- 
pendent, shifting from Vb — at fc^T » €d to Vb = V* 
at ksT <C e d - 
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FIG. 14: The shifted position of the mid-charge point V* 
versus e d , for T L /D = Tb/D — 0.01. Here A = 2.3, while 
the number of NRG states retained is equal to N s = 2300. 
Solid line: The weak-coupling formula of Eq. 1)43 ^ . Inset: 
An enlarged image of the low-e^ regime. The dashed line is 
merely a guide for the eye. 

Throughout this section we focus on a symmetric level, 
€d = —U/2, for which the position of the degeneracy 
point is pinned at Vb = for all U. The discussion of 
the asymmetric case, €d ^ is left for a forthcoming 

publication. 



VI. INTERACTING LEVEL 

Thus far, our discussion was restricted to a noninter- 
acting level. However, any realistic setup is bound to 
include an on-site repulsion on the level. If the level is 
realized by a small quantum dot, this on-site repulsion 
U will in fact exceed the charging energy Ec, due to the 
reduced geometry of the smaller dot. From the stand- 
point of the two-channel Kondo effect, the inclusion of 
a finite U raises several basic questions. Primarily, does 
the two-channel Kondo effect persist for an interacting 
level? Since an on-site repulsion couples the two spin 
orientations on the level, it is not at all clear whether 
the physical spin still acts as a passive spectator for the 
screening of the charge fluctuations inside the quantum 
box. Let us suppose that it does, what is then the effect of 
a finite U on the two-channel Kondo temperature? Can 
one still have a non-exponential Tr-? Finally, as is well 
known, an on-site repulsion U can itself introduce non- 
trivial many-body physics by forming a local magnetic 
moment on the level. Can there be a novel interplay be- 
tween single- channel screening of the magnetic moment 
on the level and two-channel overscreening of the charge 
fluctuations inside the quantum box? 

In this section we answer these questions, first at the 
qualitative level within the framework of perturbative 
scaling, and then in a detailed manner using the NRG. 



A. Weak coupling 

Similar to the case of a noninteracting level, one can 
exploit the weak-coupling limit Tl ^> , U to gain some 
analytic insight as to the effect of a finite U. Starting 
from the U — Hamiltonian of Eq. (|18|l . the effect of a 
nonzero U is to introduce the interaction term 

Hu = U Jde 1 Jde 2 Jde 3 Jde i W(e u e 2 ,e 3 ,e i ) (56) 



where 



W(ei, Ca , ea, e 4 ) = ^Jpf (ei)pf (c 2 )pf (c 3 )pf (e 4 ). (57) 

Here p^ ff (e) is the effective isospin-up DOS of Eq. I|21(l . 
while the normal orderings in Eq. H56(l account for the 
single-particle energy of the level, = —U/2 [note that 
we do not include the latter energy within the Hamilto- 
nian of Eq. {HJ]. Replacing p c ^ (e) with the symmetric 
rectangular DOS of Eq. the Hamiltonian of Eq. fT%|) 
reduces to that of Eq. (|25|l , while the interaction term of 
Eq. I|56(l simplifies to 



U 



D L rD L r-D L j-D L 

de x I de 2 / ^£3 / de 4 

D L J~D L J-D L J-D L 



(58) 
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To treat the resulting Hamiltonian, we resort to per- 
turbative scaling. Similar to the case of a noninteracting 
level, this is done in two stages. First the larger band- 
width of the box is scaled down from D to Dl = 
to arrive at an effective Hamiltonian with a single joint 
bandwidth Dl- This step is not affected by the presence 
of a nonzero U, reproducing the Hamiltonian of Eq. i|28|) 
with the extra interaction term of Eq. I|58|l . Once a sin- 
gle bandwidth is reached, one can proceed with conven- 
tional RG. Successively reducing the single bandwidth 
using poor-man's scaling, one notes the following prop- 
erties: 

• Of the two isospin Kondo couplings J±_ and J z , only 
the renormalization of J z is directly affected by the 
interaction term of Eq. 

• All interaction terms generated under the RG con- 
serve the physical spin, containing an equal number 
of creation and annihilation operators for each spin 
orientation. In particular, aside from the renormal- 
ization to Jj_ and J z , the isospin-exchange interac- 
tion retains its original two-channel form. 

• The Coulomb interaction U has a scaling dimen- 
sion of two, and is formally irrelevant. The same is 
true of all higher order interaction terms generated 
under the RG, involving four fermion operators or 
more. 

Thus, all higher order interaction terms tend to dimin- 
ish as the bandwidth is reduced, leaving only the three 
interaction terms Eq. (|28[1 : J±, J z , and V. We therefore 
conclude that the two-channel Kondo effect is robust for 
3> r# against the inclusion of a weak on-site repul- 
sion U . The latter has the effect of only weakly modifying 
the Kondo temperature Tk- 



B. Large on-site repulsion 

Another limit largely tractable by analytical means is 
that of a large repulsion on the level, U ^> Tl, Tb, e\ Vb\- 
In this limit, a stable local moment is formed on the level 
as the temperature is reduced below U/2. Therefore, 
one can resort to a generalized Schrieffer- Wolff trans- 
formation^ to eliminate the charge fluctuations on the 
level. This produces an effective low-energy Hamiltonian, 
which can be treated in turn using perturbative RG. 

To implement this strategy, we go back to the Hamil- 
tonian of Eqs. ©-© and ©, and consider the limit of 
a large repulsion on the level, U > D. Carrying out 
a Schrieffer- Wolff-type transformation, a host of spin- 
exchange ® isospin-exchange interactions are generated. 
These are conveniently expressed in terms of two sets of 
operators, acting as spin-exchange and potential scatter- 
ing from the standpoint of the physical spin: 
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(/!, v = 0,x, y, z). Here we have represented the magnetic 
moment on the level by the spin-^ operator 



(61) 



(<r being the Pauli matrices), and identified the isospin 
indices L and B with isospin-up and isospin-down labels. 
Normal ordering of the conduction-electron operators is 
with respect to the filled Fermi seas of the lead and the 
box. We further use the notation by which Sq is the unity 
operator in the isospin space, and cr is the 2x2 unity 
matrix. In terms of the above operators, the resulting 
Hamiltonian is compactly expressed as 



T~Csw 



E 



e kaC kacr C ka(T - eV B S z 



(62) 



+JoOq + J z Oq z + J±_ 
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(63) 
(64) 
(65) 
(66) 
(67) 



Physically, the spin-exchange interactions Oq and Oq z 
do not scatter electrons between the lead and the box, 
conserving thereby the isospin component S z . By con- 
trast, the spin-exchange interaction O s xx + O s yy involves 
hopping between the lead and the box. This separation of 
processes is analogous to that of intra-lead coupling and 
inter-lead coupling in the conventional two-lead Kondo 
problem. The additional O zz and O zo terms are equiva- 
lent in turn to the J z and V interactions in the Hamilto- 
nian of Eq. (J2SJ. 

Proceeding with perturbative RG, no new interaction 
terms are generated to second order as the bandwidth 
is reduced. Converting to the dimensionless couplings 
Jfi = pJfi and V v = pV v with /i = 0, z, _L and u = 0, z 
(for simplicity we assume a common density of states p 
for the lead and the box), the dimensionless couplings 
evolve according to the RG equations 



djp 
dl 



Jo + Jz 
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FIG. 15: The scaled capacitance at the degeneracy point, 
C(0, T)/cbTk /e 2 , versus T, for different values of the on-site 
repulsion U. Here T L /D = V B /D = 0.0039, A = 2.3, and the 
number of NRG states retained is equal to N 3 — 2300. With 
increasing U , a sharp shoulder develops in C(0, T) just prior 
to the onset of the characteristic two-channel logarithmic 
temperature dependence, O(0, T) ~ {e 2 /20k B T K ) ln(T fc /T). 
Concurrent with the development of the sharp shoulder, the 
onset of the ln(T) temperature dependence is pushed down in 
temperature from T/T K ~ 0.2 for U < Tl, to T/T K ~ 0.02 

for [/>r L . 



-£-J j„ 

dJ± ~ ~ l~~ 
~dT = ^ + 2^ 
dl4 _ 3 ~ 2 

dV 



dl 



= 0. 



(69) 
(70) 
(71) 
(72) 



Here I is equal to ln(D/D'), where D 1 is the running 
bandwidth and D is the bare bandwidth. 

Although Eqs. i|68|) - (|72[l have no simple analytic so- 
lution, one can read off the essential physics from the 
structure of these equations and the initial conditions 
of Eqs. H63(l - (|67|l . Primarily, the system flows toward 
strong coupling for any ratio of Yl to Tb, indicating the 
emergence of a Kondo effect for any Tl,Tb <C U. For 
Tl ~ Ts, the coupling Jj_ is the largest throughout the 
RG flow, and is the first to become of order unity. Hence 
the magnetic moment f and the isospin S are simulta- 
neously quenched. By contrast, for either Yl S> Tb or 
Tb ^> Tl, the couplings J7b and J% both grow to order 
unity at some characteristic temperature T* , while J7j_ 
and V z remain small at this temperature. In this case, 
the isospin S remains essentially free when the magnetic 
moment f is screened either by the lead electrons (for 
T L > T B ) or by the box electrons (for T B > T L ). A 
second-stage quenching of S is expected at some lower 
temperature, yet this low-temperature regime lies beyond 



FIG. 16: The two-channel Kondo temperature Tk versus U, 
for Yl/D = Tb/D = 0.0039. Here A = 2.3, while the num- 
ber of NRG states retained is equal to N 3 = 2300. With in- 
creasing U, the Kondo temperature initially grows, reaching a 
maximum for U/Tl ~ 9, and then decays. For U 3> Tl + Tb, 
when the level is in the local-moment regime, Tk can be fit- 
ted by the exponential form UbTk/D — Aexp(—BU/FL), 
where A — 0.055 and B — 0.19 (solid line). Inset: Conver- 
gence of Tk with the number of NRG states retained, for 
U/Tl = 8.9 (filled squares) and U/Fl = 16.5 (empty circles). 
For U/Tl = 8.9, T K is converged to within 2% above 2200 
states. For U /Tl = 16.5, Tk still varies by some 25% in going 
from N s = 2300 to A^ s = 2800. 



the range of validity of Eqs. H68|) - (|72|l . Below we confirm 
the two-stage quenching of f and S using the NRG. 

It should be emphasized that, regardless of the ratio 
Tl/Tb, the above analysis is insufficient for determining 
the precise nature of the low-temperature fixed point, 
whether the isospin moment is exactly screened or over- 
screened. As we show below using the NRG, the low- 
temperature fixed point is indeed that of an overscreened 
isospin moment, for all values of Tl, Tb ^ 0. 



C. General on-site repulsion 

Following the analytic analysis presented above, we 
now turn to a systematic study of all coupling regimes 
using the NRG. As emphasized in sec. IV Al it is straight- 
forward to incorporate a nonzero repulsion U within the 
NRG, as it only enters the local Hamiltonian Hq. The 
computational effort for treating a nonzero U remains 
the same as that for a noninteracting level. 

Analyzing the finite-size spectra generated by the 
NRG, wc find that the low-temperature fixed point re- 
mains that of the two-channel Kondo effect, for all values 
of U, Tl, and Tb explored. Accordingly, the capacitance 
C(0,T) diverges logarithmically with T — > for all val- 
ues of U, as demonstrated in Fig. ^| for Tl = Tb- The 
effect of a finite U in Fig. ED is most clearly seen in the 
crossover regime, prior to the onset of the characteristic 
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FIG. 17: The two-channel Kondo temperature Tk versus U, 
for Y B /D = 0.0039 and different ratios of F L to Tb- Here 
A = 2.3, while the number of NRG states retained is equal 
to N s = 2300. With increasing T L /r B > 1, the peak in T K 
as a function of U becomes shallower, until it disappears. In 
particular, no peak is left for Tl/Tb = 4. By contrast, the 
peak sharpens with decreasing Tl/Tb < 1. 



two-channel logarithmic temperature dependence of the 
capacitance. With increasing U, a sharp shoulder devel- 
ops in C(0, T) just above the onset of the ln(T) temper- 
ature dependence, which in turn is pushed down in tem- 
perature from T/T K ~ 0.2 for U < T L to T/T K ~ 0.02 
for U ^> r^. Here Tk is the two-channel Kondo tem- 
perature, extracted from a logarithmic fit to Eq. (JS2J. 
As we show below, the sharp shoulder that develops in 
C(0, T) is a signature of the simultaneous quenching of 
the isospin S and the local magnetic moment f that forms 
on the level for a large U. It is lost for Tl ^S> Tb, when 
r is quenched well ahead of the isospin S (see Fig, IT^|) . 

Figure 1161 shows the two-channel Kondo temperature 
T K versus U, for T L /D = Tb/D = 0.0039. Quite sur- 
prisingly, the Kondo temperature initially grows with in- 
creasing U, reaching a maximum for U/Tl ~ 9. This 
regime of enhanced Tk is neither covered by the weak- 
coupling analysis of sec. IVI Al nor by the large-?/ treat- 
ment of sec. lVIBI Note, however, that despite the three- 
fold enhancement of Tk as compared to the U = case, 
the logarithmic temperature dependence of the capaci- 
tance sets in at a lower temperature for U /Tl « 9 than 
for U = 0. For a large on-site repulsion, U » Tl + Tb, 
the Kondo temperature decays exponentially with U. 

Repeating the calculation of Tk versus U for different 
ratios of Tl to Tb, we find a qualitative difference be- 
tween T L > T B and T B > T L . As seen in Fig. El the 
peak in Tk as a function of U becomes shallower with in- 
creasing Tl /Tb > 1) until it disappears. For Tl/Tb = 4, 
for example, there is no peak left. The slope of Tk versus 
U at U — also becomes shallower as Tl/Tb is increased, 
in agreement with the perturbative scaling analysis of 
sec. IVI Al The latter predicts a weak U dependence of 



FIG. 18: The capacitance (unsaturated curves) versus the 
magnetic susceptibility of the level (saturated curves), for 
Tl/D — Tb/D — 0.0039 and different values of the on-site 
repulsion U. Here A = 2.3, while the number of NRG states 
retained is equal to N 3 = 2300. Both the capacitance and 
the magnetic susceptibility are scaled with the two-channel 
Kondo temperature Tk, extracted from the slope of the ln(T) 
diverging term in the capacitance. Explicitly, ItsTa/D is 
equal to 0.0025, 0.0047, and 0.0019 for U/D = 0.01, 0.034, and 
0.07, respectively. For U/D = 0.01, the level is on the boarder 
line between the mixed-valent and local-moment regimes. As 
U is increased, a stable local moment forms on the level. In 
this regime, the onset of the ln(T) temperature dependence 
of the capacitance and the saturation of the magnetic suscep- 
tibility take place roughly at the same crossover temperature. 



the Kondo temperature for U, Tb -C Tl . 

In contrast, the peak in Tk versus U becomes sharper 
and higher as Tl /Tb < 1 is decreased. For Tl/Tb = 1/4, 
Tk actually exceeds the conduction-electron bandwidth 
as one approaches the peak position. In this range of U, 
there is a clear separation (over three orders of magni- 
tude) between the low-temperature scale Tk, extracted 
from the slope of the ln(T) component of C(0,T), and 
the thermodynamic crossover scale To, below which the 
ln(T) temperature dependence sets in. In fact, it becomes 
increasingly difficult to extract a meaningful Tk from the 
slope of the ln(T) diverging term as one approaches the 
peak position, possibly signaling the breakdown of the 
two-channel Kondo effect at some critical U. We em- 
phasize, however, that the low-temperature fixed point 
revealed by the NRG remains that of the two-channel 
Kondo effect, for all values of U explored. 

Up to this point, we focused our discussion on the 
overscreening of the isospin S, as probed by the capaci- 
tance. However, for a large U, a local magnetic moment is 
formed on the level. To clarify the interplay between the 
screening of the magnetic moment and the overscreening 
of the isospin S, we compare in Figs.ll8landll9lthe capac- 
itance (i.e., the isospin susceptibility) with the magnetic 
susceptibility of the level. To this end, we augment the 
Hamiltonian of Eqs. © (JSJ) and J3J| with the local mag- 
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FIG. 19: The capacitance (dashed lines) versus the mag- 
netic susceptibility of the level (solid lines), for U/D — 0.3, 
Tl/D — 0.0157, and different ratios of Tl to Tb- Here 
A = 2.3, while the number of NRG states retained is equal to 
N s = 2300. The capacitance is scaled with the two-channel 
Kondo temperature Tk, extracted from the slope of its ln(T) 
diverging term. The magnetic susceptibility is scaled with 
its zero-temperature value, which defines yet another Kondo 
scale fcsTich = (/isg) 2 /4x(0). For Fl = Tb, the magnetic 
susceptibility saturates at the same temperature range where 
the ln(T) temperature dependence of the capacitance sets in. 
As the ratio Tl/Tb is increased, the impurity moment is 
quenched well ahead of the isospin S. Accordingly, a large 
separation builds up between the associated one-channel and 
two-channel Kondo scales fcsTi c h and \zbTk (marked by the 
full and dashed arrows, respectively). 



netic field 



7~Lrnag — — ^B9jHt z , 



(73) 



where t z is the z component of the spin operator defined 
in Eq. (|61|l . The magnetic susceptibility of the level is 
given in turn by the derivative 



X(T) = fi B g. 



d(r z ) 
dH 



(74) 



evaluated at zero field. 

Figure El shows our results for Tl/D — Tb/ D = 
0.0039. In the local-moment regime, the magnetic sus- 
ceptibility saturates at low temperatures, much in the 
same way as it does in the conventional one-channel 
Kondo effect. In accordance with the large-Z7 analysis 
of sec. IVIBI the saturation of x(T) occurs at the same 
temperature range as the onset of the ln(T) tempera- 
ture dependence of C(0, T), confirming the simultane- 
ous quenching of the spin r and the isospin S. Indeed, 
the zero-temperature susceptibility x(0) is of the order of 
(HBgj) 2 /Tk, indicating that the same Kondo scale Tk 
underlies both the screening of f and the overscreening 
of S. 



As the ratio Tl/Tb is increased, see Fig. El two 
distinct Kondo scales emerge, one associated with the 
single-channel screening of r, and the other associated 
with the two-channel overscreening of S. As before, the 
two-channel Kondo temperature Tk is extracted from a 
logarithmic fit of the capacitance to Eq. (|54ll . while the 
one-channel Kondo temperature is defined from the zero- 
temperature magnetic susceptibility: 



knT 



(PBfj) 2 

4x(0) 



(75) 



For T L = r B (we fix T L at T L /D = 0.0157 and vary 
T B in Fig. the two scales T K and T lch are essentially 
the same, differing by a factor of three in favor of Tk- 
Upon decreasing by a factor of four, i.e., Tb — Tl/^, 
the order of scales is reversed, and Ti c h becomes twenty- 
five- fold larger than Tk- Upon further reducing to 
Tb = r^/16, the scale Ti c h exceeds Tk by four orders of 
magnitude. Hence the impurity moment r is quenched 
well ahead of the isospin S for Tl ^S> Tb, in accordance 
with the analysis of sec. I VI Bl In this limit one can safely 
speak of two successive Kondo effects: first the impu- 
rity spin undergoes one-channel screening by the lead 
electrons, followed by two-channel overscreening of the 
charge fluctuations inside the box. A similar qualitative 
picture is recovered for Tb Tl, except that the impu- 
rity spin is screened by the box electrons rather than the 
lead electrons. 

Interestingly, the sharp shoulder that characterizes the 
capacitance for Tl = Tb is lost for Tl S> Tb, and one 
recovers a capacitance line shape that closely resembles 
the U — case. In particular, the onset of the ln(T) tem- 
perature dependence of the capacitance is pushed back 
to T/T K ~ 0.2 in Fig.El The sharp shoulder that devel- 
ops in C(0,T) for Tl = Tb <C U is therefore a distinct 
signature of the simultaneous quenching of S and f in 
this case. 



VII. DISCUSSION AND CONCLUSIONS 

The charging of a quantum box, weakly connected to a 
lead by a single-mode point contact, is one of the leading 
scenarios for the realization of the two-channel Kondo 
effect. 4 The main obstacle hampering the observation of 
this effect in semiconductor quantum dots stems from the 
exponential smallness of the Kondo temperature Tk- In 
order for a fully developed two-channel Kondo effect to 
emerge, ksTK must significantly exceed the level spac- 
ing inside the box. However, with Tk being exponentially 
smaller than the charging energy Eq, it is practically im- 
possible to realize a measurable Kondo scale that exceeds 
the level spacing in present-day semiconducting devices & 

In this paper we proposed a possible remedy to this 
problem, by considering a setting in which tunneling be- 
tween the lead and the box takes place via a single reso- 
nant level. The basic idea is to exploit the strong energy 
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dependence of the transmission coefficient through the 
impurity to achieve nearly perfect transmission at the 
Fermi energy, but only small transmission away from the 
Fermi energy. In this manner, the large transmission at 
the Fermi energy enhances Tk , while the small transmis- 
sion away from the Fermi energy insures the emergence of 
a sharp Coulomb-blockade staircased This should be con- 
trasted with the case of an energy-independent transmis- 
sion coefficient, where the Coulomb staircase is washed 
out for perfect transmission^ 

As seen in Figs. El and El the Kondo temperature 
is indeed dramatically enhanced when the impurity is 
tuned close to perfect transmission at the Fermi energy. 
Specifically, Tk varies approximately linearly with the 
level broadening T l + T b , and is many- fold larger than 
the exponentially small Kondo scale obtained for a tun- 
nel barrier with comparable tunneling matrix elements. 
We emphasize that this enhancement of Tk is not re- 
stricted to perfect alignment of the level with the Fermi 
energy. Rather, it extends also to a level off resonance 
with the Fermi energy, provided that the single-particle 
transmission coefficient is large at the Fermi energy. This 
point is demonstrated in Fig. El where IzbTk / (Tl + T b) 
roughly equals 0.2, even though is as large as Tl = Tb- 
Note that the corresponding transmission coefficient at 
the Fermi level is indeed large, being equal to T = 0.8. 

For T < 1, the Kondo temperature is again small. 
However, we find a qualitative difference between <C 
and Tb ^> T^. For Tb <C Fl, the quantum box is 
only weakly coupled to the level, and thus to the lead. 
Accordingly, the Kondo temperature takes the exponen- 



tial form Tk oc exp 



-Tr 2 /2y/T , similar to the case of a 

weak tunnel barrier with a transmission coefficient equal 
to TA This result was derived both for a level at reso- 
nance with the Fermi energy, Eq. I|30|l , and for a level off 
resonance with the Fermi energy, see end of sec. IIVBI 

By contrast, the Kondo scale deviates from the ex- 
ponential dependence on 1/ VT, for Tb ^> Ft,. For 
| erf | ^> Tb Tl, for example, Tk depends in a power- 
law fashion on T^, as specified in Eq. I|44|l . For a level 
at resonance with the Fermi energy, the ratio Tb/T^ de- 
fines the crossover from weak (Tb/Tl -C 1) to strong 
(Tb/Tl 3> 1) coupling, which emphasizes the inequiva- 
lent roles of Tl and in the two-channel Kondo effect 
that develops. 

Any practical realization of a resonant level will nec- 
essarily involve an on-site Coulomb repulsion U on the 
level. Such an interaction couples the two spin orien- 
tations, questioning the development of the two-channel 
Kondo effect. In this paper we focused on a symmet- 
ric level, U + 2e<j = 0, which is easier to study numeri- 
cally since both the degeneracy point and the mid-charge 
point are pinned at Vb = 0. For an asymmetric level, 
U + 2cd 7^ 0, the degeneracy point (assuming it exists) 
is shifted away from Vb = 0, and is no longer identified 
with the mid-charge point. For U = 0, this was explicitly 
shown to be the case in Figs. 1121 and 1131 

As seen in Fig. El the two-channel Kondo effect is 



robust against the inclusion of an on-site repulsion, for 
a symmetric level. For Tl = Tb, the Kondo scale ex- 
tracted from the slope of the hi(T) diverging term in 
the capacitance is actually enhanced by a moderately 
large repulsion U (see Fig. fT5)l . although the onset of 
the logarithmic temperature dependence is pushed down 
to lower temperature. For T l < Tb and intermediate 
U , there are clearly two distinct energy scales govern- 
ing the low-temperature capacitance: a thermodynamic 
crossover scale &bT , below which the ln(T) temperature 
dependence sets in, and a low-temperature scale IzbTk, 
which fixes the slope of the m(T) diverging term. This 
differs from the conventional two-channel Kondo Hamil- 
tonian, where To ~ Q.2Tk and Tk are essentially the 
same. 

For a large on-site repulsion, U ^> Tl + Tb, the Kondo 
temperature is again exponentially small, this time in 
U/{T L +T B ), see Figs. El and El For T L > r s , one 
can qualitatively understand this result along the fol- 
lowing lines. As seen in Fig. El f° r Tl ^> Tb there 
is a sequence of Kondo effects: first the local moment 
on the level undergoes single-channel screening by the 
lead electrons, followed by two-channel overscreening of 
the charge fluctuations inside the box. Below the single- 
channel Kondo temperature Ti c h oc exp [— ttU/STl], the 
spectral function of the d electrons consists of two high- 
energy resonances at ±U/2, and a narrow Abrikosov- 
Suhl resonance of width Ti c h and height 1/ttTl, which 
is pinned at the chemical potential^ It is the latter 
peak that serves as the effective lead density of states, 
pf{e) of Eq. (SSJ, available for screening the charge fluc- 
tuations inside the box. Substituting T lc h in for in 
Eq. (|29|l , the pre-exponential factor in Eq. is reduced 
by a factor of \fT\ c \y/Y l, yielding a Kondo temperature 
which is exponentially small in both U/Tl and ^JTl/T b- 
Obviously, the above picture overlooks the possible rel- 
evance of higher order interaction terms generated upon 
the screening of the d local moment. 

In semiconducting devices, one can realize a tunable 
level using an ultrasmall quantum dot, whose charging 
energy U is bound to exceed that of the quantum box. 
Since the bandwidth D in the effective Hamiltonian of 
Eqs. ©-© is °f the order of Ec, this dictates the hier- 
archy U > D ^> Tl,Tb- Hence, not much can be gained 
from a symmetric level, as the associated Kondo temper- 
ature Tk is exponentially small in U/{Tl + Tb) 1. 
Instead, one would like to tune the level to the mixed- 
valent regime, where Tk is expected to be of the order of 
the level broadening. 

As emphasized above, treatment of such an asymmet- 
ric level is complicated by the need to accurately locate 
the position of the degeneracy point, which is no longer 
pinned at Vb = 0, and does not coincide with the mid- 
charge point. In fact, one cannot entirely rule out the 
possibility that the two-channel Kondo effect is unstable 
against particle-hole asymmetry for an interacting level, 
as recently implied by Le Hur and Simon. 33 Using per- 
turbative RG and an analogy to a related model of two 
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capacitively coupled quantum dots^a these authors ar- 
gued that the two-channel Kondo is lost in the local- 
moment regime, when a stable local moment is formed 
on the levels Although no explicit distinction was made 
between a particle-hole symmetric and an asymmetric 
level, the analysis of Le Hur and Simon implicitly as- 
sumed an asymmetric level, by taking the coupling V to 
be nonzero (see Ref . l33f) . For V = 0, there are qualitative 
changes to the RG flow of Le Hur and Simon. Indeed, 
Figs. 1 151 1 191 unambiguously establish the emergence of 
the two-channel Kondo effect for 26^ + U = 0, including 
in the local- moment regime U Tl,Tb- 

For weak coupling, Tl ^> Tg, |e<j|, U, one can analyti- 
cally see that particle-hole asymmetry does not play any 
role in the emergence of the two-channel Kondo effect. 
To this end, we extend the analysis of sec. IVI Al to an 
asymmetric level, 2e^ + U ^ 0. The effect of a nonzero 
Ae = Ed + U/2 is to supplement the Hamiltonian terms of 
Eqs. (|18fl and (|56l) with an additional potential-scattering 
term for the a\ L(J degrees of freedom: 

Hps = Ae^2 Jde 1 Jde2v(e 1 , e 2 ) :a\ lLa a^La- (76) 

with 

v(e 1 ,e 2 ) = ^/pf(e l )pf(e 2 ). (77) 

Replacing p°^(e) with the symmetric rectangular DOS of 
Eq. 1)24(1 , the potential-scattering term of Eq. 1)76(1 reads 

Hps = -f^y]/ de i de 2 :al La a e2La : , (78) 



which is just the sum of the two Hamiltonian terms of 
Eq. (ED) with D m -► D L and V+ = V- = Ae/(27rr L ). 
Aside from renormalizing the voltage Vp , the addition of 
such a Hamiltonian term does not affect the low-energy 
physics. The system continues to show the two-channel 
Kondo effect, albeit at a shifted position of the degen- 
eracy point. This differs from the conclusion of Le Hur 
and Simon for the local-moment regime. 

Away from weak coupling, a full-scale numerical effort 
is required to resolve the effect of particle-hole asymme- 
try on the emergence of the two-channel Kondo effect. 
Such a study is currently under way. Our preliminary 
results for the mixed-valent regime (setting td equal to 
zero) support the conclusion that the two-channel Kondo 
effect is robust against particle-hole asymmetry for an 
interacting level. Moreover, there are indications that 
Tk remains significantly enhanced in the mixed-valent 
regime also for U several times larger than D ~ Ec, 
as is the case in realistic semiconducting devices. A de- 
tailed analysis of this physically relevant regime will be 
presented in a forthcoming publication. 
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